Dynamic positive column in long-gap barrier discharges. 
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A simple analytical model of the barrier discharge in a long gap between opposing plane electrodes 
is developed. It is shown that the plasma density becomes uniform over large part of the gap in the 
course of the discharge development, so that one can speak of a formation of a dynamic positive 
column. The column completely controls the dynamics of the barrier discharge and determines such 
characteristics as the discharge current, discharge duration, light output, etc. Using the proposed 
model, all discharge parameters can be easily evaluated. 

PACS numbers: 



Oh! 



CZ3 
O ■ 

43 ■ 



<n : 

o . 
Tt- ■ 
o ■ 
in ; 
o . 

o : 

>>: 
43 : 

Or 



X 



I. INTRODUCTION 

Dielectric barrier discharges are typically generated in 
gas gaps between two electrodes covered with dielectric 
layers. Sinusoidal or square-wave voltage with frequen- 
cies ranging from few kHz to hundreds of kHz are used to 
generate the discharges. The product of the gas pressure 
and the distance between the dielectric surfaces pD can 
vary quite significantly. In the present paper we consider 
barrier micro-discharges such as used in Plasma Display 
Panels Q and excimer lamps where pD ps 5 — 50. 

The dynamics of a barrier discharge during one pulse of 
the applied square-wave voltage was considered in detail 
analytically in our previous work [j| , where it was shown 
that at high overvoltage the discharge develops into an 
ionization wave, whose velocity is determined primarily 
by the charge production rate in the cathode fall region. 
This wave moves from the anode toward the cathode, 
resulting in contraction of the cathode fall region and in- 
crease of the electric field within this region. Upon reach- 
ing the cathode, the ionization wave can either quickly 
disappear (when the capacitance of the dielectric layer is 
small) or transform into a quasi-stationary DC cathode 
fall (when the capacitance is large). The main assump- 
tion of that work was that the resistance of the plasma 
trail, created by the ionization wave, can be neglected. 

In the present paper we consider the opposite case, 
when the barrier discharge dynamics is strongly influ- 
enced by the plasma trail which electrically connects the 
anode with the cathode fall region. As will be shown be- 
low, in the case of a long gap the plasma trail becomes 
uniform over the large part of the gap, virtually forming 
a dynamic positive column. The uniformity of the col- 
umn enables us to consider it simply as a variable resistor 
through which the cathode fall (CF) charges the dielec- 
tric layer capacitor (see Fig.0. When the capacitance is 
large, one can introduce further simplifications into the 
model. In this case the CF is, in essence, quasi-stationary 
and its V-I characteristics can be approximated by those 




FIG. 1: Equivalent circuit of a long-gap barrier discharge. 

of the DC cathode fall Q. 

In our work we neglect the volume and near-wall re- 
combination of charged particles, which if included into 
consideration would only increase the resistance of the 
positive column and enhance its role in the barrier dis- 
charge dynamics. 

Although dynamic characteristics of barrier discharges 
have been already extensively studied via computer sim- 
ulations (see Ref. || and references therein), it remained 
unclear, for example, what processes control the duration 
and amplitude of the current pulse. The proposed model 
answers this question and can be used to estimate var- 
ious parameters of long-gap barrier discharges between 
opposing plane electrodes. This type of discharges (as 
well as |a 00) ma y prove to be a viable alternative to 
near-surface discharges 0, 0, 0, favored nowadays 
in Plasma Display Panel industry. 



II. QUALITATIVE CONSIDERATION AND 
BASIC EQUATIONS 

We will consider discharges in a gas gap between op- 
posing plane electrodes covered by dielectric layers of 
thickness d and dielectric constant e. We assume that 
the gap length L g is much larger than the length of the 
normal cathode fall, 



Lg > £„ 



(1) 
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It is also assumed that the dielectric layer capacitance 
is large, i.e. the effective thickness of dielectric layers is 



2 



small, 



2d/e < L n 



(2) 



For example, in the case of 10%-Xe/90%-Ne mixture at 
pressure 500 Torr and secondary emission coefficients 
iNe = 0.5 and jxe = 0.005, the normal cathode fall 
length Lnorm ~ 6.5/zm, the gap length under considera- 
tion is from several hundred microns to one millimeter, 
and the effective thickness of dielectric layers is less than 
or about several microns. 

Although discharges in gaps with such a long distance 
between the anode and cathode can be initiated in sev- 
eral different ways (by, for example, using a set of aux- 
iliary electrodes), we assume that the Townsend break- 
down takes place . 

Under these assumptions (long gap, large capacitance, 
Townsend mechanism of the breakdown) , the general pic- 
ture of the discharge dynamics can be described as fol- 
lows: 

1. Upon application of a sufficiently high voltage across 
the gap (C/ ap pi > £/br, where Ubr is the breakdown volt- 
age), the Townsend criterion is satisfied and the positive 
charge starts to build up in the gap. While the amount of 
total charge is small (and thus the electric field is undis- 
turbed), it grows exponentially with time. 

2. At some moment it reaches a critical value (~ 
€oVa. P pi/Lg) and causes considerable distortion of the 
electric field. The field vanishes at the anode and from 
this moment on there coexist two different regions in the 
gap: the region filled with plasma (plasma trail) where 
the electric field is relatively small, and the region of the 
gap adjacent to the cathode where the electric field is 
strong and the electron density is negligible. With time, 
the plasma trail expands toward the cathode while the 
CF region contracts. At this stage, the process of the 
plasma trail uniformization begins. 

3. When the length of the CF region becomes com- 
parable to inorm, it transforms into a quasi-stationary 
DC cathode fall and the discharge current sharply in- 
creases and, after reaching its maximum, decreases. The 
plasma density in the already uniformized plasma trail 
(positive column) grows while the electric field there falls 
(see Fig.|5J). These processes are accompanied by the in- 
tensive deposition of power from the external source into 
the discharge. At the end of this stage, the main portion 
of the external voltage drops across the dielectric layer 
capacitor, which leads to quenching of the discharge. 

4. In the afterglow, the number of ions and electrons 
in the gap gradually decreases through the dissociative 
recombination of electrons and molecular ions. Charged 
particles are also pulled out of the plasma toward the 
dielectric surfaces by the residual (and ambipolar) elec- 
tric fields. Note that although the discharge current is 
small during the discharge decay, the voltage transferred 
to the capacitor is about the normal CF voltage V noim 
(when 2d/e <C L noim ). 

The key feature of long-gap barrier discharges is that 
the plasma trail is uniform during the stage three, when 
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FIG. 2: Electric field E and electron density n e in the gap 
at time moments ti < ti < t$ when the discharge current is 
strong: I(ti) = 0.5I ma x, = /max, and I(t 3 ) = 0.5/ max , 

where J max is the maximum discharge current. 



the discharge current is strong. The intrinsic mechanism 
responsible for the uniformization can be qualitatively 
explained as follows: Due to the plasma quasi-neutrality, 
which quickly establishes in the Maxwellian time (see 
Sec. II Vj) . the particle current is constant throughout the 
plasma trail. Therefore the electric field is higher in re- 
gions where plasma density is lower. The higher electric 
field results in higher ionization rate, which in turn leads 
to faster growth of the plasma in these regions and, even- 
tually, to a leveling off of the plasma density. 

To find the plasma density in the positive column, one 
can use the continuity equation for electrons 



dn dnjj, e E 



Of 



dx 



r(E) n, 



(3) 



where n is the electron (plasma) density, /i e is the elec- 
tron mobility, E is the electric field in the column, and 
r(E) is the ionization rate. For simplicity, we omitted 
the diffusion part of the electron flux in Eq. © . Neglect- 
ing the motion of ions (ion mobility is several hundred 
times lower than the electron mobility) and using the 
constancy of the particle current, one can immediately 
arrive at the following relationship between the plasma 
density and the electric field in the positive column, 



en/2 e E = I(t)/S, 



(4) 



where e is the elementary charge, I(t) is the discharge 
current, and S is the cross-sectional area of the discharge. 
Since the electron mobility is almost constant at small 
electric fields, from Eq. it follows that 



SnE pc + n pc 5E 



0, 



(5) 
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where Sn and SE are perturbations and n pc and E pc are 
the average values of the plasma density and electric field 
in the plasma column. Using Eqs. 0-0; one can readily 
obtain 



9n E 



dt 



dSn 



(E pc ) - E, 



dr(E pc ) 



pc 



dE, 



pc 



(6) 
(7) 



While the obvious equation I© describes growth of the 
plasma density in the positive column, Eq. {7\ shows fast 
decay of the density perturbations (since under positive 
column conditions E pc r'(E pc ) > r(E pc )). Let us em- 
phasize that, as is seen from our derivation, the mech- 
anism of the positive column uniformization has an en- 
tirely dynamic nature and manifests itself regardless of 
the diffusion processes. It is interesting to note that the 
uniformization takes place even in a more general case 
[l3| when the discharge can not be treated as quasi-onc- 
dimensional. 

In this paper we give a simple analytic description of 
only the third stage of the discharge development, when 
discharge current is strong. During this stage, the posi- 
tive column occupies almost the entire gap, and as seen 
from Eq. © its resistance R w L g / (en pc n e S) is governed 
by 



dR , s 
= -r{E pc ) R, 



where 



E n 



IR 



(8) 



(9) 



The quasi-stationarity of the cathode fall allows us to use 
the V-I characteristic of the corresponding DC cathode 
fall, 



V cf = Vcf(I). 



(10) 



The applied voltage is distributed between the elements 
of the circuit in Fig. ^ according to the second Kirchoff 
law, 



IR + V CF (I) + ^ = V apph 



(11) 



where C = eoeS/2d is the capacitance of the dielectric 
layer capacitor, and the capacitor charge q changes with 
time as 



dq = 
dt 



(12) 



The equations JBJ-lJEJ) form the foundation of our 
model of barrier discharges with long positive columns 
and a quasi-stationary cathode fall. They must be com- 
plemented with the following initial conditions 



R(t - 
q(t 



-00 J = 00, 

-oo) = 0. 



(13) 
(14) 



Note, that at not very high current densities, the DC 
CF voltage changes insignificantly (see Fig. |3J| . The case 
when the Vqf is taken constant and equal to V nomTl is 
considered in the next section. 
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FIG. 3: V-I characteristics of the DC cathode fall obtained 
in the fluid approximation for the gas mixture and secondary 
emission coefficients mentioned at the beginning of Sec. [H] 
Dotted curve corresponds to the unstable subnormal cathode 
fall. 



III. CONSTANT CF VOLTAGE 
APPROXIMATION 



When V CF = K, 



const, it is convenient to use 



R and E pc as independent variables. Replacing the IR 
with E pc L g and after that taking the time derivative of 
Eq. I|ll|) , one can obtain the equation for the electric field 
in the positive column, 

dE _ E 
~dt = ~~RC' 



(15) 



where the subscript "pc" is omitted for brevity. Eqs. © 
and (|15|) must be solved with the following initial condi- 
tions 



E(t 



R(t 

Eq 



■ — oo) 

Vappl — 



= 00, 
V CF 



(16) 
(17) 



Condition (|16|l indicates that at the beginning, there is no 
plasma in the discharge gap, and condition (|17|l directly 
follows from condition l|14|l and the Kirchoff 's law Qlip. 

We would like to point out that, based on what was 
said in the previous section, Eqs. JSJ and (|T3|l do not 
describe the initial stages of the discharge development, 
when the positive column is being formed, neither do they 
describe the final stage, when the discharge quenches and 
the cathode fall dies. 

Now let us turn our attention to the analysis of the 
system of Eqs. JHJ and l|15|) . Dividing left and right parts 
of equation JSJ) by the corresponding parts of equation 
(|15f) and multiplying both parts by C we obtain 



d{RC) r{E) 



dE 



E 



{RCf 



(18) 
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FIG. 4: Typical dependence of I and E on time. 



Integrating this equation with initial conditions (| 1 fcil) and 
(fPTjl gives 



RC 



where 



rEQ 

f(E,E ) = \ 

JE 



Kg) 

E 



dE. 



(19) 



(20) 



Substitution of (|19f) into l|15f) leads to a differential equa- 
tion for the electric field, 



dE_ 
~dt 



-Ef(E,E ), 



(21) 



which should be solved with the initial condition Q17[l. 
From this equation, one can see that the time evolution of 
the electric field in the neutral column depends only on its 
initial value and the ionization rate, and does not depend 
explicitly on the electron mobility. Because the right 
part of equation (|21l) is always negative, the electric field 
monotonically decreases with time from Eq to zero. The 
discharge current in our model is uniquely determined by 
the electric field, 



ELg 

R 



= LgCE f(E,E ). 



(22) 



It goes from zero (at E = Eq), trough the maximum (at 
E = E*), to zero again (at E = 0), see Fig.0J The elec- 
tric field E* at the current maximum can be determined 
from the extremum condition dI(E)/dE — 0, which, with 
the use of Eq. (|22|l . can be written as 



r(E*)=f(E*,E ). 



(23) 



The typical shape of the function r(E)/E (which is pro- 
portional to the the ionization coefficient a(E)) as well as 
the graphical solution of equation i)23[l is given in Fig. [5] 
From this figure, one can see, in particular, that always 
Eq/2 < E* < E . If r(E)/E grows faster than the linear 
function, then E /^ < E* < Eq. 

Using Eqs. (1221) and 1|23[1 . the expression for maximum 
current can be written as 



CL g E* r(E* 



(24) 




FIG. 5: Graphical solution of equation 1231 . Area of the 
dashed rectangle on the left is equal to r(E t ), while the area 
of the dashed region under the r(E)/E curve on the right 
is equal to f(E,,Eo). These areas are equal when E* is the 
electric field at the current maximum. 



This equation shows that in the present model the max- 
imum current is proportional to the dielectric layer ca- 
pacitance. Dependence on the gap length L g is more 
complicated, since E* also depends on L g (through Eq). 

The time duration r cur of the current pulse can be de- 
fined as the ratio of transferred charge Qtr = C (Kippi — 
Vcf) = CLgEo to the maximum current, 



En 



E„ r(E*) 



(25) 



It is useful to note that when r(E)/E grows faster than 
the linear function, then 



> 



5.2 



3\/3 
r(Eo) ~ r(E )' 



(26) 



For a steeper dependence of the ionization coefficient on 
the electric field, the numerical coefficient in Eq. I|26l) is 
even greater. 

From (|19|l immediately follow the expressions for the 
final resistance of the positive column, 



R(t 



1 



oo 



Cf(0,E o )' 

and for the final plasma density in the column, 

CL 



n(t = oo) 



r(0,E ) 



(27) 



(28) 



Using the expression for the ionization coefficient a = 
r/[i e E, the ratio of the total number of electrons (ions) 
created in the positive column N — nL g S to the number 
of transferred electrons (ions) Qtr/e can be written in the 
form 



eN 



= (a) L a 



(29) 



where 



E 



E0 



a{E) dE. 
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Note that (a) < 0.5a(£'o) if a{E) grows faster than 
the linear function. If, in addition to that, V app i < 
Vb r + V no tm, then there is a practically useful estimate 
for the number of charged particles created in the posi- 
tive column, 

eN < ^a(E )L g Q tI < X - In (l + -)CV ap pi. (30) 

This expression gives reasonable estimates far beyond the 
applicability of our model. Note that since the electric 
field is small in the positive column, in the case of noble 
gas mixtures, mainly the ions of species with the lower 
ionization potential are created there. 

Similar formulae can be obtained for the number of 
excited atoms (molecules) created in the positive column 
during the discharge, 

^7T^ = {a C xc)L g , (31) 

1 f EO r cxc (E) 
(cw) = — / - dE, 

where r cxc is the excitation rate. Again, in the case of no- 
ble gas mixtures, mostly the species with lower excitation 
energies get excited. 

Using Eqs. and the energy deposited into the 
positive column can be estimated as 

W pc = ^C{V appl -V CF ) 2 , (32) 

while the total energy which goes from the external 
source into the discharge is 

Wtot = \CV^ X . (33) 

Note, that the same amount of energy, CV^ pl /2, is stored 
in the capacitor and is recovered during the next dis- 
charge pulse. 

The difference between Eqs. Q33JI and (|32J) gives the 
amount of energy that goes (in the considered approxi- 
mation) into the cathode fall, 

Wcf = \c[Vl pl - (K pp] - V CF ) 2 } . (34) 

IV. DISCUSSION 

To verify the validity of the approximation Vcf = 
Unorm, one should substitute I max from Eq. (|2~4"|) into the 
V-I characteristics of the DC cathode fall (|10|l and make 
certain that 

^CiK-fmax) ~ U norm . (35) 

If the obtained CF voltage is noticeably higher than the 
normal CF voltage, one should analyze the full system of 
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FIG. 6: Renormalized discharge current for L g — 800/im 
and Vkppi = 390V: (1) Vcf = Unorm approximation; (2) 
fluid simulation for 2d/e = 1/im; (3), (4) fluid simulation 
and Vcf = Vcf (I) approximation for 2d/e — O.lfim; (5) 
fluid simulation for 2d/e = 10/^m. Time is measured from the 
moment of maximum of the currents. 
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FIG. 7: The same as in Fig. |S1 but for L g — 800/im and 
U appl = 450U. 
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FIG. 8: The same as in Fig. HJ but for L g — 400/im and 
U appl = 320U. 

Eqs. ||8J|-m), i.e. use the more accurate Vcf = Vcf [I) 
approximation instead of the V = U norm approximation. 
In Figs. 161^1 are shown results obtained from our model 
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and from fluid simulations under different discharge con- 
ditions. A striking feature of all these figures is that the 
amplitude of the renormalized current I/C changes by 
only a factor of 2 — 3 when the capacitance C changes by 
two orders of magnitude (compare curves 2,3 and 5). A 
prediction of our model that I/C does not depend on C 
(see Eq. (|23J) is qualitatively correct even when the cri- 
terion © of the model applicability is not satisfied. For 
curve 5, the effective dielectric thickness 2d/e = lO/im is 
greater than the normal CF length L norm ~ 6.5/im, so 
that the CF is not quasi-stationary but dynamic! 

Understandably the Vcf — Vnorm approximation over- 
estimates the amplitude of the discharge current, since 
Vnorm is the minimal voltage on the cathode fall; this 
overestimation is more pronounced for shorter gaps. 

Another interesting feature of the long-gap discharge 
dynamics is that the duration of the current pulse in- 
creases with increase of the gap length (compare Fig. |S| 
with Figs. El and Q, which is not at all obvious since the 
applied voltage is greater in the case of the longer gap. 
Qualitatively this discharge property is attributed to the 
increase in inertia of the positive column with its length. 

Note also that in the case of the large capacitance (ab- 
normal quasi-stationary CF, see curve 3), the discharge 
current grows quickly and falls slowly. In the opposite 
case of the small capacitance (dynamic CF, see curve 5), 
the current grows somewhat slower than it falls. 

Similar to DC discharges, the positive column in bar- 
rier discharges is more efficient source of light (radiation) 
than the cathode fall. In the positive column, where the 
electric field is weak, almost all the deposited energy goes 
into excitation of atoms, while in the cathode fall signif- 
icant portion of the energy is spent on the ion heating, 
and, because the electric field is strong, a considerable 
part of the rest goes into ionization of atoms rather than 
excitation. It is obvious that the case when the cathode 
fall is very abnormal (i.e. when the C is very large) is 
not optimal for the light production. The question of the 
optimal capacitance of dielectric layers - i.e. the optimal 
relationship between 2d/e and L n orm (which determines 
to what degree the CF is dynamic) - should be dealt with 
in the context of concrete discharge parameters. 

At the end of this section, let us make several re- 
marks regarding the quasi-neutrality of the plasma in 



the positive column and possible generalizations of our 
model. The quasi-neutrality results from the fact that 
Maxwcllian relaxation time is much less than the time of 
plasma growth in the column (which is of the order of 
r cur ). Using Eq. l(2"5|) . one can estimate 

eo eo'S' 2d 

•7" = i 

e\i e n CL g f(0,E o ) eL g f(0,E o )' 

Since f (0, Eo) = r(0, EL) + f{E*,E ) > r(J5») ~ r"^ the 
quasi-neutrality condition takes the form 



7"cur 

Our model allows us to include plasma losses due to am- 
bipolar diffusion to the walls surrounding the discharge 
volume and due to recombination by subtracting from 
the r.h.s. of Eq. © the terms tT^ n pc (raf is the char- 
acteristic ambipolar diffusion time) and k r rip C (k r is the 
recombination coefficient). Note, however, that the side 
walls can strongly influence the breakdown stage of the 
discharge 01 an d can even lead to the formation of near- 
dielectric-surface discharge rather than the volume one. 



V. CONCLUSION 

The dynamics of a barrier discharge between oppos- 
ing electrodes separated by a long gap which is filled 
with a mixture of noble gases has been considered. It 
has been shown that during the discharge development, 
there forms a region between the anode and the cathode 
fall similar to a DC positive column, where the plasma 
density and electric field are uniform but, in contrast to a 
DC case, change in time. The mechanism of the positive 
column uniformization has an entirely dynamic nature. 
It is caused by the ionization processes and has nothing 
to do with the diffusion. 

A simple discharge model allowed us to obtain useful 
estimates for the plasma density, number of excitations 
and ionizations, and power consumption; it also helped 
to capture non-trivial trends in the dependence of the 
discharge current on discharge parameters. 
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